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I. INTRODUCTION 

This report sunnnarizes the work done thus far xmder Contract NAS5-20974 
concerning development of the three-axis stabilized navigation model for 
purpose of measuring winds from the ATS-6 image data. 

Over the past year our work with the ATS-6 attitude telemetry data 
and the associated image data has shown that the telemetry data is not of 
sufficient quality for the purposes of measuring winds accurate to a few 
meters per second. However, a methodology has been developed to account 
for the relative attitude changes between successive ATS-6 images which 
allows reasonable high-quality wind sets to be produced. The method consists 
of measuring the displacements of the right and left infrared earth edges 
between successive ATS-6 images as a function of scan line; from these 
measurements the attitude changes can be deduced and used to correct the 
apparent cloud-displacement measurements. The preprocessing time and effort 
for obtaining these earth-edge measurements is comparable to that of obtaining 
landmark measurements. 

The wind data sets generated from ATS-6 using the "earth-edge" metliodology 
mentioned above were compared with those derived from the SMS-1 images (and model) 
covering the same time period. Quantitative comparisons for low-level 
trade cumuli were made at interpolated uniformly spaced grid points and for 
selected individual comparison clouds. There were sixteen grid points at 
which vector comparisons could be made between the ATS-6 and SMS-1 derived 
winds. The differences between each satellite-derived wind vector component 
at each of these points were all less than 2 meters per second with a 
root-mean-square difference of less than 1 meter per second for both vector 
components. Wind direction differences did not exceed 13®. Tor the nine 
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selected individual comparison clouds, the root-mean-square differences for 
the U and V components were 1.0 and 1.2 meters per second with a maximum 
wind direction difference of 15°. 

^ t 

The major task remaining imder this contract is the implementation of 
the ATS-6 model (including the new earth-edge methodology) to AOIPS system 
at GODDARD; this work is currently in progress. 
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II. ATS-6/SMS-1 CLOUD VELOCITY COMPARISON MEASUREl^TS 

« 

This sectioa discusses the ATS-6 wind sets produced using the methodology 
discussed briefly in the INTRODUCTION and in more detail in SEC. Ill and 
con 5 )ares them to the SMS-1 derived winds. 

The background on the study area and the procedure used to obtain 
both ATS-6 and SMS-1 wind sets using McIDAS are given in the 31 June, 1976, 
NAS5-20974 progress report. The SMS-1 wind sets and all associated figures 
and measurements discussed in that report are reproduced here unchanged 
except for figure and table captions. 

FIGURES II. 1. and II. 2. show the ATS-6 and SMS— 1 low-level trade cumuli 

wind fields and FIGURES II. 3. and II. 4. their associated wind fields 

interpolated to 2® latitude/longitude grid points. TABLE II. 1. summarizes 

the low-level wind-vector comparisons for small clouds which could be 

identified on both the SMS-1 and ATS-6 images, and TABLE II. 2. summarizes 

the low-level wind-vector comparisons for the interpolated grid-point 

values shown in FIGS. 11^.3. and II. 4. Note that the root-mean-square values 

in the component differences (AU , AV ) are all about 1 m/sec for both 

rms rms 

the small comparison clouds and the interpolated grid-point values. The 
slight bias of the ATS-6 V-components to values smaller than the corresponding 
SMS-1 values for the small-cloud comparisons (TABLE II. 2.) is curious because 
these same clouds are scattered throughout the clouds in the study area 
from which the interpolated grid-values are obtained — and the interpolated 
ATS-6 V-component values are slightly higher on the average than their SMS 

counterparts. This small-cloud bias most likely results from the small 
number of clouds measured. 

FIGURES II. 5. and II. 6. are plots of the high-level ATS-6 and SMS-1 
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'winds. ■ Quantitative comparisons were not made because of the small number 
of measurements available. 

FIGURE II. 7. is a plot of the wind-vector differences of the successive 

pairs of ATS-6 images and serves as a quality check on the vectors. These 

vectors are related directly to those given in FIGURE II. 1., where in this 

figure, the U-component, for example, was obtained by averaging U and U 

12 23 

where U^^ is the U-component measured for a cloud between images T^ (164223Z) 

and T (170658Z) and U between images T and T (173134Z). In FIG. II. 7. 

2 23 2 3 

however, the U-component is (U - U )/2 and therefore is a measure of the 

12 23 

distance from the mean of these two values. The same discussion applies 
to the V-component. Thus, FIG. II. 7. is an indication of the consistency 
and accuracy of the ATS-6 wind-measuring process. Ideally these vector 
differences would be zero for non-accelerating clouds (effectively true for 
this study area), and if not zero, at least small and pointing in random 
directions. This latter condition is generally satisfied as implied by 
FIG. II. 7. TABLE II. 3. presents the vector differences for the two ATS-6 
image pairs for the same small-cloud measurements given in TABLE II. 1. and 
also makes comparisons to the SMS-1 values. The standard deviations given 
in TABLE III. 3. show that there is slightly less variation between the ATS-6 
wind-vector component measurements and ^^^(''^23^ there is between 

^12^^12^ and the corresponding averaged SMS components Furthermore, 

of the ATS-6 vector-component differences are less than 0.5 m/sec. 


the averages 
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DAY 74195 TIME170658 900-900 MB WINDS 


FIGURE II. 1. LOW-LEVEL WINDS FROM ATS-6 computed by using earth-edge displacement 
measurements to account for attitude changes. Length of vectors are proportional 
to velocity. Reference vector shown in lower left-hand corner. 



DRY 74195 TIME170000 900-900 MB WINDS 


FIGURE II. 2. LOW-LEVEL 
shown in FIG. II. 1. SMS- 
ATS-6 due to the former 
lower left-hand comer. 


WINDS FROM SMS-1 covering same time period as ATS-6 wind field 
•1 wind field is of smaller geographical extent than than of 
satellite's higher resolution. Reference vector is shown in 
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FIGURE II. 3, LOW-LEVEL ATS-6 WINDS INTEEIOLATED AT 2® GRID POINTS from the data shown in 
FIG. II. 1, Vector components at each grid point for which both ATS-6 and SMS-1 (FIG. II. 4) 
non-zero vectors are available are tabulated and compared in TABLE II. 2. Reference vector 
is shown in lower left-hand corner. 
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No. 


SMS 

SMS 


ATS- 

■A 

DIFFERENCE 




LAT 

LON 

U(m/s) 

V(m/s) 

U(m/s) 

V(m/ s) 

AU 

AV 

IV%I 

A^ 

1 

22.99 

-66.75 

- 2.82 

4.25 

- 1.37 

4.14 

1.45 

-0.11 

1.45 

15 

2 

21.75 

-67.48 

- 4.14 

5.27 

- 4.42 

3.94 

-0.28 

-1.33 

1.36 

-10 

3 

21.12 

-67.76 

- 3.99 

5.15 

- 3.28 

2.61 

0.71 

-2.54 

2.64 

-14 

4 

20.10 

-63.46 

- 7.45 

2.51 

- 7.65 

1.98 

-0.20 

-0.53 

0.57 

- 4 

5 

18.21 

-62.81 

- 9.14 

2.31 

- 7.87 

2.16 

1.27 

-0.15 

1.28 

1 

6 

16.10 

-74.88 

-10.95 

3.76 

- 9.48 

2.64 

1.47 

-1.12 

1.85 

- 3 

7 

16.05 

-67.13 

-11.22 

1.36 

-11.08 

0.15 

0.14 

-1.21 

1.22 

- 6 

8 

15.97 

-63. 84 

-11.00 

-2.28 

-12.25 

-3.76 

-1.25 

-1.48 

1.94 

- 5 

9 

15.17 

-64.34 

-11.16 

-1.44 

-11.79 

-1.47 

. -0.63 

-0.03 

0.63 

0 








AU=0.30 
S(AU)=0.'98 
AU =0.97 

^=-0494 
S(AV)= 0.82 
AV = 1.22 




9S 
^ 2 

!?3 ^ 

|S_ 


TABLE II. 1. ATS-6/SMS-1 LOW-LEVEL WIND-VECTOR COlsff^ISONS FOR SMALL CLOUDS which could be identified 
on both the SMS-1 and ATS-6 images. The averages X and standard deviations S _^bo^h in m/s) are given 
for the U- and V-component differences between the two sets of wind vectors; |v^-Vgj is the magnitude 
of their vector differences and A(j) the differences in wind direction. 
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SMS-1 


ATS-6 


DIFFERENCE 


LAT 

LON 

U (m/s) 

V (m/s) 

U (m/s) 

V (m/s) 

AU 

AV 

IV''sl 

A(j) 

23 

-70 

- 3.09 

5.18 

- 2.31 

li.37 

.78 

- ,81 

1.12 

3° 

23 . 

-68 

- 2.67 

3.89 

- 2.28 

3.96 

.39 

.07 

.40 

5° 

23 

-66 

- 2.46 

3.72 

- 1.59 

3.79 

1 

.87 

.07 

.87 ' 

11 ° 

23 

-62 

- 4.30 

3.24 

- 4.26 

3.98 

.04 

.74 

.74 

6 ° 

21 

-68 

- 4.16 

5.00 

- 3.65 

3.15 

.51 

-1.85 

1.92 

- 9° 

21 

-64 

- 5.93 

2.13 

- 7.03 

j 1 

1.61 

- 1.10 

- .52 

1.22 

- 7° 

21 

-62 

- 7.14 

2.22 

- 8.88 

1.85 

-1.74 
1 1 

- .37 

1.78 

- 6 ° 

19 

-64 

- 7.21 

1.83 

- 7|.21 

2.06 

1 

.00 

.23 

.23 

2 ° 

19 

-62 

- 8.56 

1.15 

- 9.05 

2.05 

! 1 1 

- .49 

.90 

1.02 

5° 

17 

-78 

- 8.03 

- .25 

- 6.71 

0.37 

1 

1.32 

.62 

1.46 

5° 

17 

-76 

- 8.06 

1.79 

- 7.80 

3.44 

.26 

1.65 

1.67 

11 ° 

17 

-74 ■ 

- 9.15 

2.08 

- 9.03 

2.69 

.12 

.61 

.62 

4° 

17 

-72 

'''-•9.24 

1.08 

- 8.27 

i 

2.75 

.97 

1.67 

1.93 

12 ° 

15 

-78 

-11.60 

- .06 

i 

-10.81 

1.15 

.79 

1.21 

1.45 

6° 

15 

-76 

, -10.49 

1.70 

-10.47 

1. 31 

.02 

- .39 

.39 

- 2 ° 

15 , 

’ 1 , 

.-64^ 

—10 . 13 

-1.99 

-11.61 

-1.13 

-1;48 

.86 

1.71 

6 ° 


AU=0.09 AV=0.29 

• S(AU)=0.88 S(AV)=0.94 

AU • =0.86 AV =0.95 
rms rms 

TABLE 11,2. ATS-6/SMS-1 LOW-LEVEL WIND-VECTOR COMPARISONS FOR INTERPOLATED GRID-POINT VALUES. 
Vector components at each grid point for which both ATS— ^ (FIG. 11.3.) and SMS— 1 (FIG. II, 4.) 
non— zero vectors are available are given. The averages X and standard deviations S (both in m/s) 
are eiven for the U- and V-component differences. 
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DAY 74195 TIME170658 200-300 MB WINDS 

FIGURE II. 6. HIGH-LEVEL SMS-1 WINDS derived from same imagery as SMS-1 low-level 







DRY 74195 TIME170658 900-900 MB WINDS 


FIGURE II. 7. ATS-6 WIND-VECTOR DIFFERENCES. Vector components are (D 12 - U23>/2 
and (Vi 2 - V23)/2 where U . , , V.. are the velocity components measured from the 
ith - jth image pair. TheSe values are tabulated in TABLE II. 3. and compared to 
SMS^l values. Double circles in figure indicate vectors of zero magnitude. 





No 


Jo. 

ATS- 

± 


SMS 


ATS- 

-A 


SMS 

1 


Ui2(m/s) 

U23(m/s) 

AU, 

A 

^S 

AD 

S 

Vi2<m/s) 

V23(m/s) 


^S 


1 

- 2.07 

- 0.67 

-0.70 

- 2.82 

0.75 

5.29 

3.00 

1.15 

4.25 

1.04 

2 

- 3.84 

- 5.00 

0.58 

- 4.14 

0.30 

4.85 

3.04 

0.91 

5.27 

-0.42 

3 

- 4.03 

- 2.54 

-0.75 

- 3.99 

-0.04 

4.22 

1.01 

1.61 

5.15 

-0.93 

4 

- 6.71 

- 8.60 

0.95 

- 7.45 

0.74 

1.51 

2.45 

-0.47 

2.51 

-1.00 

5 

- 8.16 

- 7.60 

-0.28 

- 9.14 

0.98 

2.05 

2.30 

-0.13 

2.31 

-0.26 

6 

- 9.19 

9.80 

0.31 

-10.95 

1.76 

3.06 

2.24 

0.41 

3.76 

-0.70 

7 

-12.51 

- 9.66 

-1.43 

-11.22 

-1.29 

0.04 

0.28 

-0.12 

1.36 

-1.32 

8 

-12.52 

-12.00 

-0.26 

-11.00 

-1.52 

-4.59 

-2.95 

-0.82 

-2.2g 

-2.30 

9 

-11.13 

-12.46 

0.67 

-11.16 

0.03 

-0.84 

-2.11 

0.64 

-1.44 

0.60 




X— 0.10 
S=0.79 


X=0.19 

S=1.06 



X=0.35 

S=0.80 


X=-0.59 

S=1.00 


TABLE II. 3. WIND- VECTOR DIFFERENCES OF SUCCESSIVE PAIRS OF ATS-6 IMAGES and comparison to SMS-1 values. 
These data are for the same nine small clouds given in TABLE II. 1. U^2 U— component measured from 

ATS-6 images Ti and T2, U23 for images T2 and T3. Au» = (U12 - U23>/2 is the difference from their mean 
value and AU_ = Ui 2 “ difference of U12 from the mean SM^-1 measured U— component Ug for the same 

cloud. The same deflnifions apply to the V-components. Averages X and standard deviations S given under 
the difference columns are in meters per second. 
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“III. METHODOLOGY TO ACCOUNT FOR ATTITUDE CHANGES WHEN MEASURING AND 
COMPUTING WINDS 

1. Method and Resiilts 

The method for removing the attitude changes between successive ATS-6 
images by using earth-edge measurements is described below. 

For simplicity let us assume that wind measurements are to be made 
between only two successively scanned ATS-6 images designated Tj and T 2 . 

The following procedure is followed: ^ 

(i) The two infrared images are loaded on successive TV frames on McIDAS 
with the same start-coordinates so chosen that the right earth-edges can be viewed 
and the scan lines cover the same scan-line range' cooresponding to the study area. 
Using an image-matching technique, the displacement of the T 2 ~edge relative 

to the T^-edge is measured along the scan lines in increments of approximately 
10 TV-scan lines (5 ATS-6 scan lines) at a time until measurements have been 
made along the entire available earth edge. These data are recorded and 
curve fitted (see FIG. III.l.), It is important to point out that the image- 
matching technique used measures displacements only along scan lines and 
not perpendicular to it. 

(ii) The same procedure as given in (i) is performed for the left 

I 

earth edge corresponding to the same scan line range. These measurements 
result in data such as shown in FIG. III. 2. 

(iii) It is shown in Appendix A that from these right and left earth- 
edge displacement curves it is possible to coD 5 )ute the satellite displacements 
(AL(L), AE(L)) of the T 2 - SSP (Sub-Satellite Point) relative to the Ti - SSP 
as a function of mirror-scan line number L. (The resulting AL(L) curve and 
AE(L) curve derived from the curves presented in FIGS. III.l. and III. 2. are 
shown in FIGS. III. 3. and III. 4. respectively). The AL curve is in effect 
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RIGHT EDGE IR T1 - T2 



MIRR0R SCRN 


FIGURE III.l. RIGHT EARTH-EDGE DISPLACEMENT MEASUREMENTS as a function of mirror- 
scan number. The crosses show the actual measured displacements (measured by an 
image-matching technique) in resolution elements for a particular mirror-scan line 
of the T 2 image relative to the Ti image of the right earth edge. The infrared 
digital imagery was used to obtain these measurements. The measured data were fit 
with a 15th degree polynomial curve shown in the figure. 
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ELEMENT DISPLRCEMENT 


LEFT EDGE IR T1 - T2 



MIRROR SCAN 


FIGURE III. 2. LEFT EARTH-EDGE DISPLACEMENT IIEASTJREMENTS . These measurements are 
complementary to those shown in FIG. III.I. 




LINE SHIFT 


i' 




FIGURE III. 3. RELATIVE LINE SHIFT OF THE T 2 SSP RELATIVE to the Ti-SSP as a function 
of mirror scan number. This curve was derived from the curves sho^vrn in FIGS. III.l. 
and III. 2. by the methods discussed in the text and Appendix A. (SSP = Sub-Satellite 
Point; Tj corresponds to ATS-6 image DAY 195 1974, 164223Z and T 2 , 170658Z same day) 
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ELEMENT DISPLRCEMENT FR0M 15 DEGREE CURVE 



FIGURE III. 4. RELATIVE ELEMENT SHIFT OF THE T 2 -SSP RELATIVE to the Ti-SSP (solid 
curve) as a function of mirror scan number- This curve is the complement to the 
one in FIG. III. 3. Dotted curve is pitch-angle difference between images Ti and T 2 
obtained from the attitude telemetry data plotted in FIGURE III- 5- and clearly 
shows that these data are not able to account for the actual observed changes 
between AIS-6 images. 
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the roll-angle difference curve between Tj and T 2 as a function of line 
ntmiber and the AE curve the pitch-angle difference curve between these two 
Images. 

(iv) The AL, AE curves are converted to look-up tables as a function 
of line number. Thereafter, whenever a cloud displacement is measured 
between the T^-and T 2 -images for the purposes of wind velocity measurements, 
the AL, AE values corresponding to the average scan line number of the cloud 
are "looked-up" and these values are subtracted off the cloud-displacement 
coordinates obtained by the conventional image-matching technique. The 
resulting modified vectors are then run through the ATS-6 model as before 
to convert satellite image coordinates to earth-coordinates and hence wind 
vectors. That this simple procedure does in fact remove the roll and pitch 
angle changes correctly is justified in APPENDIX B. 

The method discussed above has the following advantages : 

(i) attitude telemetry data are not needed; 

(ii) because the method is basically a null-detection scheme, the 
"earth-edge" need not be defined; 

(iii) distortions in the scan system do not affect the nature of the 
wind fields or the magnitude of the wind vectors; 

(iv) landmark measurements are required for only one image in a study 
sequence (this statement is clarified in Appendices B and C) . 

The reason the method works is basically because: 

(i) the ATS-6 attitude control system during the image acquisition time 
was good enough such that attitude changes were generally small and therefore 
small-angle app.roximations can be used in the ATS-6 navigation model with a 
high degree of accuracy; 
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(il) attitude changes are generally slow compared to the time required 
to scan a line (period for an attitude oscillation cycle is on the order 
of 100 mirror-scan lines ) ; 

(iii) distortions in the mirror-scan nonlinearities are small and 
slowly varying (+2 resolution element change over 200 resolution elements 
along a scan line - see Appendix G) ; 

(iv) mirror-scan nonlinearities repeat in magnitude and phase from 
line to line and change very slowly from image to image; 

(v) yaw-angle changes are small (see Appendix E) ; 

(vi) eccentricity of ATS-6 orbit is very small (see Appendix F) . 

2.,^ Attitude Telemetry Data 

FIGURES III. 5-7 show the pitch, roll and yaw-angle differences between 
the Tj to T 2 images used to obtain the wind-measurements reported in this 
report. These attitude data were obtained from the NASA/ GODDARD magnetic 
tapes supplied to us and represent the latest processing procedure applied 
to the ATS-6 telemetry data. Our June, 1976, progress report showed plots 
of the attitude data derived from a different processing technique which 
contained "saw-teeth" in the yaw and pitch data; it is evident that this 
artifact does not appear in the present data. That report also presented 
evidence that the large pitch and roll changes around Mirror Scan 200 were 
indeed "present" in the associated image data but that the smaller attitude 
angle changes seemed to have little relationship with the studied image 
data. The same statement is true of the present attitude data. The Tj - T 2 
pitch-angle differences in FIG. III. 5. corresponding to the Study Area 
(from which the wind measurements were obtained) have been plotted in 
FIG. III. 4. (dotted curve). Recalling that the Element Shift curve (solid 
line) in this figure is essentially the correct -pitch-angle-difference curve. 
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FIGURE III. 5. PITCH-ANGLE DIFFERENCES (DEGREES) BETT^EEN ATS-6 IMAGES Ti (164223Z) 
and T 2 (17065 8Z) attitude telemetry data from latest NASA/ GODDARD tapes supplied 
to us. Data corresponding to Study Area mirror-scan range was plotted in FIG. III. 4. 
(A first-order pitch difference results in an element difference between successive 
ATS-6 images; a 0.0083° pitch change results in a 1-element shift) 
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T1-T2 R0LL DIFFERENCES. NEN DRTH 



FIGURE III. 6. ROLL-ANGLE DIFFERENCES (DEGREES) for same image pairs discussed in 
FIG. III. 5. 


ORIGINAL PAGE IS 
OP POOR QUALny 







YflW DIFFERENCES 


III-IO 



FIGURE III, 7. YAW-ANGLE DIFFERENCES (DEGREES) for same image pairs discussed in 
FIG. III. 5. 









it is apparent that attitude telemetry data contains no useful information 
for the purposes of accurate wind determination. If the dotted instead 
of the solid line had been used to correct for attitude changes when < 
computing the reported winds, errors of +5 m/sec would have resulted. 

Comparing the roll-difference plot (FIG. III. 6.)- to the Line Shift curve 
(FIG. III. 3.), the opposite situation exists in the sense that the roll- 
difference change in the Study Area interval suggested by the attitude data 
(about 1.5 m/sec) was not in fact observed to be the case (the Line Shift 
curve is essentially plot over this mirror-scan range) . 

3. Error Analysis 

There are two potential-error sources related to attitude changes which 
could reduce wind measurement accuracies : 

(i) large yaw-angle changes (discussed in Appendix E) ; 

(ii) large amplitude and frequency roll-angle changes occurring in 
the vicinity of +^00 ATS-6 scan lines about the SSP-scan line (discussed 
in Appendix B) , 

For the ATS-6 image data which is so afflicted, there is nothing that 
we know about that can be done to correct for these affects. i 
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APPENDIX A. ANALYSIS OF IIETHOD TOICH DETERMINES ATS-6 SSP IMAGE COORDINATE 
DISPLACEMENTS BETlffiEN SUCCESSIVE IMAGES RESULTING FROM ATTITUDE 
CHANGES 

As shown in Section III, the ATS-6 attitude telemetry data time series was 
shown to be inadequate for the purpose of accurate wind measurements. This 
appendix presents an analysis of the technique used to calculate the ATS-6 
Subsatellite Point (SSP) image coordinate changes between successive images 
(shown in FIGS. III. 3. and III. 4.) resulting from attitude changes during the 
image-scan time. Appendix B shows how these measurements are used to account 
for the attitude changes in the process of computing accurate cloud displacements. 

1. Method 

Let 'Tj , T 2 designate two successive ATS-6 data images where Tj^ is the 

"earlier" of the two images and (L , E ) (Line scan and Element numbers) 

c c 

are its SSP image coordinates (FIG. A. 1.) determined from the ATS-6 navigation 
model. Let (L,E) be the T^^ image coordinates for a point on the right 
earth edge and AE , AE the measured displacements along line L of the 

XC Li 

f 

Right and Left T^-earth edges relative to the T^^ earth edges. (These 

measurements are obtained in practice on McIDAS using the infrared ATS-6 

data images and an image-matching technique which is constrained to measure 

displacements of the earth edge only along a scan line. It is worth 

emphasizing that the image-matching method indeed measures displacements 

of the geometrical earth edge and not features near it - such as clouds.* 

This is not surprising since the greatest contrast is between earth and 

space - not within features near the earth's edge.) 

The object, then, is to compute the displacement coordinates (AL,AE) 

of the T 2 -SSP relative to the T^-SSP at line L using the measured values 

of AE ,AE . 

K L 

*See added note at the end of this appendix. 




FIGURE A.l. Earth-edge Displacement Measurement 
Geometry 



A-3 


- Pj = 2 V(P^ - P^,)2 - (Pj - p^,)2 


= 2 VTp - p ■) 2 Vp _ p \2 
^3 O'' ^5 0’^ 


= 2 V(L - L^)2 + (e - E^)2 - [L - (L^ + L)]2 , and (3) 

= 2(E - E^) . (4) 

Substitution of (3) and (4) into (2) and rearranging yields 

6 =, Vx2 + y2 - (Y - AL)^ - X , where (5) 

6 = (AEj^ - AE^)/2 
X = E - E 

c 

Y = L - L . 

c 

Since 6 is obtained from the AE_,AE^ measurements and (L,E), (L ,E ) are 

K L c c 

known, AL can be solved for (5) to yield 

AL = Y + \^AD , where (6) 

RAD=y2-2X5-62, ' 


The choice of sign in (6) is determined by substituting the expression for 
6 in Eq. (5) in RAD; the result is 


AL = Y + V(Y - AL)^ = Y + , Y > AL 

= Y - VrS) , Y < AL , 


(7) 


The peculiar condition implied by (7) that AL must be known before AL can 
be calculated is really not a problem since generally AL is small (+ 5 lines) 
and [y| is usually » |Al|. For most cases then the conditions are 

AL = Y + VrAD , Y < 0 , (8) 


•RAD=Y^-2X6-62 . 
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For the case where [y| approachs (AL) in value there are other problems; 
these are discussed at the end of this appendix. 

Eqijations (1) and (8) , then, define the displacement in image coordinates 
of the T 2 -SSP relative to the Tj^-SSP for line L. By repeating this entire 
process for other scan lines the SSP-shifts (AL,AE) as function of T^-line 
position are obtained. In practice the right and left earth-edge displace- 
ments are first measured over the entire range of scan lines of interest. 
Curves are then separately fit to the right and left edges measurements (see 
FIGS. III.l. and III. 2.) resulting in a AE„ vs L and a AE, vs L curve over 
the scan-line range of interest. The curve values themselve are then used 
in (1) and (8) to compute the AL,AE curves shown in FIGS. III. 3. and III. 4. 

2. Error Analysis 

The magnitude of the errors associated with the displacement coordinates 
(ALjAE), denoted by (A(AL), A(AE)), depend on the accuracy to which the 
horizontal earth-edge displacements can be measured and the location on the 
image (i.e. polar, equatorial) where (AL,AE) are being computed. 

For regions near the poles, the slope of the earth edge approaches zero 
and the "edge" covers an increasing number of elements for each scan line. 

The horizontal image matching accuracy does not noticeably degrade until 
the slope is below 45°, nominally corresponding to a latitude of 50°. 

However in this region of the image, the location inaccuracies resulting 
from foreshortening of features on the earth is severe enough that wind 
measurements are generally not made at these higher latitudes. Thus the 
error in AE(A(AE)) associated with the study regions of interest (between 
+50° latitude) is the error associated with the difference in the earth- 
edge displacement curves shown in FIGURES III.l. and III. 2, From (1) and (5) 
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it is seen that A(iE) = A6 and this average error as seen in these figures 
and also suggested by the accuracy in the resulting wind measurements is 
estimated to be about .1 pixel; thus, A(AE) = A6 = .1 pixel. (9) 

The error A (AL), on the other hand, increases in size for lines near 
the SSP (nominally the equator) . 

Letting 

= (L - L )^ + (E - E )^ = (radius of earth-image)^ = (1044 pixels)^ 
c c 

b^ = (Y - AL)2; Y = L - L 

c 

A6 = measxirement error, 
then from (8) 

A(AL) = AL(6 + A6) - AL(6) 

= + [VRAD(6 + A6) - VEAD(6)] for b $ 0 

^ + [Vb^ - 2A6(a^ - b2)l/2 _ b] . (10) 

The region in which equation (8) for the computation of AL breaks 
down is when the value of L (and hence b) is such that the expression of the 
radical in (8) is zero, i.e. 

Region for which (8) is invalid occurs when; 

|b| = |L - - ALj \/2aA6 = 14 pixels , (11) 

where the assumption b^ « a^ was made (generally jAL| does not exceed 
5 scan lines, thus as L approachs L^, b^ gets veiy small compared to a^). 

For b^ 5 10 (2aA6) , (10 reduces to 

A(AL) = + Va^/b^ - 1 , b $ 0 . (12) 

Values of A(AL) computed for various values of b using (11) appear 


in the table below. 
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Approximate number of 
lines from SSP-line 

|bl 


50 

100 

200 

300 

400 

500 

600 


TABLE A.l. Line- Shift Errors 


Mean-maximum error 
(lines) 

|a(al)|‘'‘ 


2.1 

1.0 

0.5 

0.3 

0.2 

0.2 

0.1 


as a Function of Line Number 


for an earth-edge displacement error of 0.1 pixel. 


'visible image lines (2 visible lines equals 1 ATS-6 mirror-sweep 
line) . 


r 


This table and TABLE D.l. indicate that in the range of + 200 lines 
(+ 100 ATS-6 scan lines) about the SSP-line that wind-measurement errors 
would exceed 1 m/sec for a 3- image study sequence if nothing were done to 
reduce the A(AL) errors. Fortunately, the AL vs L curve (FIG. III. 3.) 
is a slowly varying curve whose amplitude fluctuations tend to be small; 
thus the AL vs L curve can be Interpolated across the SSP region from 
the curves derived on either side of it. 


*Note added from p. A-1. That the image-matching technique measured earth- 
edge and not cloud displacements on the earth edge was checked by two 
methods; 1) earth-edge measurements were made on the earth's limb using 
the visible imagery and compared to those measurements made from the 
associated infrared imagery covering the same portion of the earth edge; 

2) the infrared digital Images were remapped so that brightness values 
above and below threshold values determined from the earth-edge values 
were set equal to a constant leaving only the earth edge values unchanged. 
(The unchanged earth edge image was about 3 pixels wide or 15 km at the 
equator) . Earth edge measurements were made and compared to the measurements 
made from the unmapped data. In both tests the resulting comparison curves 
were identical. 


OEIGINAL PAGE ^ 
OF POOR QUALITY 
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APPENDIX B. ANALYSIS OE ALGORITffil WHICH ACCOUNTS FOR RELATIVE ATTITUDE 
CHANGES 'IN SUCCESSIVE ATS-6 DATA IMAGES USING EARTH-EDGE 
SHIFT MEASUREMENTS 


This appendix provides a semi-rigorous analysis of the algorithm used 
to account for the relative attitude changes in a sequence of ATS-6 data 
images using the earth— edge measurements derived by the technique discussed 
in the main portion of the report and APPENDIX A. In effect, this appendix 
shows why the algorithm "works"; the following appendices discuss in more 
detail the limitations of and the errors associated with this method. 

Consider a sequence of n ATS-6 images designated by Tj^, T^, ... T^, 
from which the displacements in earth coordinates of a feature (cloud) 
between successive images are to be determined. Let 

(L.,E.) = image coordinates (Line, Element) of the cloud feature 
X X 

for the T^^^ image, 

= f^„(L,,E.) = unit pointing vector in Picture Frame coordinates 
PFx PF X X 

derived from (h^,E^) using the ATS=6 scan-camera geometry, 

= unit Local Vertical pointing vector associated with 

R, = 3x3 rotation matrix which transforms without error r^„. 

X LVx 

into fri-r.- at the time the feature was scanned, i.e. 

PFx 


^PFi \ ^LVi * 


( 1 ) 


1. A Two-Image Sequence 

Consider now the first two images , T^ in the sequence whose relation- 
ship between LV and PF coordinates are given by 


^PFl ^ \vi 

A A 

^PF2 ^ ^2 ^LV2 * 


( 3 ) 
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The relationship of the cloud feature in LV coordinates for image T^ 
can be related to as follows : 




(4) 


where 4r^^^ is the displacement of the cloud in LV coordinates from to 

T and results from two effects: 

2 

(i) apparent change in the earth’s orientation over the time 
interval caused by a non-zero angle of inclination in 

the satellite’s orbit (the eccentricity of the ATS-6 orbit 
was so small (- lO"*^) that no significant change in angular 
size occurs - see Appendix F) ; 


(ii) motion of the cloud relative to the earth’s surface. 


Now the orientation of the earth relative to the LV frame depends only on 
the orbit and the earth’s position relative to celestial coordinates; thus 
if and were accurately known, the displacement of the cloud over 

the time interval associated with these two vectors could be calculated 
with an accuracy limited only by the equations describing the dynamical 
relationship between the satellite orbit and the earth’s position relative 
to celestial coordinates. It is assumed that for ATS-6 that this transformation 
produces negligible error. 

Returning to equations (2) and (3) , the rotation matrices R^ , R^ differ 

from each other slightly because of an attitude change of the PF frame 

relative to the LV frame between T and T ; they are related to each other 

by an infinitesimal rotation ^ ^ the unit matrix and 

an antisymmetic matrix whose elements are the small-angle differences between 

the PF axes at times T, and T . Thus 

1 2 
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\ - (I + Ej,) Sj ; 


substituting (4) and (5) in (3) and expanding 


"PF2 = ^"lVI ^ ^12^"LV1 ^2^""LV^^^12> 


( 5 ) 


( 6 ) 


The first term on the RHS of (6) is the last term is negligible 

C-io'^ times smaller than the first order terms) , the second term is the change 
in PF coordinates due to an attitude change and the third term is the change 
^n PF coordinates due to the two changes in the LV frame discussed above. Using 
(4) to combine the 'first and third terms /'equation (6) then can be rewritten as 


’^PF2 ®"1^LV2 ^12^PFl * 


(7) 


Comparison of (7) with (2) implies that if and E.^ were known, 

and therefore the cloud displacement - could be computed with 

a high degree of accuracy. 

A good approximation to is obtained from the earth-edge displacement 

measurements which provide the displacement in image coordinates AL^^(L), 

t 

AE,„(L) of relative to T, as a function of line number L. The line shift 
12 2 1 

AL ^2 is proportional to a small change in roll ^hout the PF-X axis and 

the element shift is proportional to small changes in pitch ^^*^2 the 

PF-Y axis, i.e. 

‘•hi ‘ 

where p is the angular size of a pixel = 1.45 X 10”^ radians. 

^ T 

Letting r „ = [X, Y, Z] where T is the transpose, the explicit form of 

Fi* 1 

E^^ and the last term in (7) is ORIGINAL PAGE IS 

OF POOR QUALITY 

These displacement measurements are obtained by the method discussed in 
Sec. III. and Appendix A; curves showing the displacements AL,AE as a 
'funCjtion of L are shown for one image pair in FIGS, III. 3. and III. 4. 



I 


^12^PF1 


0 0 API lx 


0 0 AR Y 


-AP -AR 0 J LZ. 


= [APZ, ARZ, - (APX + ARY)] , - (9) 

where it is tanderstood that AR, AP refer to the to changes in attitude 
at the times the cloud was scanned. 

Note contains no non-zero yaw angle element (small rotation about 
the axis passing through the image SSP) since there is no way that the 
earth-edge displacement measurement technique can provide this angle; however, 
the analysis given in Appendix E shows that the effect on the accuracy of 
wind measurements (displacements) for yaw changes is generally small. Analysis 
of landmark measurements and the ATS-6 wind sets shown in this report, 
indicate that for the data images studied thus far, the yaw changes are 
negligible. 

With E ^2 and hence the last term in (7) determined, equation (7), can 
be rewritten as 

" ^PF2 “ ^12^PF1 " ^1^LV2 ’ 

where interpreted as the cloud pointing vector with the attitude 

changes removed. 

2. Computing ^ Practice 

Cloud displacement measurements and computation of ^^pp 2 accomplished 
as follows: 


(1) the image coordinates (L^,E^) of the cloud are recorded. 



B-5 


(2) the total displacement of the cloud in image coordinates 

T T 

(AL , AE ) from to are measured on McIDAS using an 
image matching technique, 

(3) from previously determined earth-edge displacement parameters, 

the displacement coordinates AL , AE due to relative 

12 12 

T T 

attitude changes are computed and subtracted from AL , AE . 
These differences are added to L^, and the associated 
PF vector computed. That this vector is a good approximation 
of r' _ is shown below. 


By definition. 


^PFl " 


( 11 ) 




( 12 ) 


Expanding £^,^(1^ + - AL^^(L^), E^ + AE^ - AE^^(^ 2 ^^ ^^1 

T 

E^ + AE ) keeping first order terms, we have 


From the geometry of the ATS-6 scan-camera, 

.T 

-SinA, COSACOSO, 

,T 


Tpy = [-cosXsinS, -sinX, cosXcosS] 


(14) 


= [X, Y, Z]- 


X = (L^ - L)p = mirror step angle, 

6 = (E - E)p = mirror sweep angle, 
c 

(L ,E ) = image center line and element value, 
c c 

p = angular size of a line or element. 


(13) 


where 
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Taking the partials of r „ in (14) with respect to L,E, and retaining 
first order terms with the partials evaluated at (L^,E^): 


3L ^PF '^^12 3E ’^PF ^^12 

= [0, PZ, -PYJ^&L + [PZ, 0, -PXl'^AE , 
12 12 ’ 


(15) 


siibstituting (8a) and (8b) into the above ei^ression and adding the two 
vectors in (15) yields 


[ZAP , ZAR , - (XAP + YAR )]^ , 

*■ 12 * 12 * ^ 12 12 ■“ ’ 

which is equivalent to (9); therefore (15) is the same as E^^r^y^. Since 
the first term in (13) is have 

r_(L +AL^ - AL (L^), E + Ae''^ - AE (L^)) = 

PF 1 12 ^ 1 12 2^ 


(16) 


r “ E r = r’ bv 
PF2 12 PF PF2 ^ 


comparison with (10). Thus the method discussed above correctly yields 


PF2* 


3. Transforming PF to LV coordinates 

With the attitude changes between T^ and T^ removed, we have 

r — R r 
PFi 1 LVl 


■PF2 


= R r. 


1 LV2 


(see (10)) 


(17) 

(18) 


Let R be the transformation from LV to PF coordinates derived from 
T^ landmark measurements in a least squares sense assuming a constant 
attitude over the time interval the T^ landmarks were scanned, i.e. R 
is a constant matrix. R then will differ slightly from R^ due to small 
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attitude changes over the scan time and, in addition, will differ from 
(assttmed to be the error-free transformation)! due- to mirror scan 
nonlinearities. Similarly, as above, we can relate R to R^^ by an infinitesimal 
transformation 

R"= R^ (X ~ e) , (19) 


where e is the infinitesimal antisymmetric error matrix which is a function 

of line and element. It is important to realize that the elements of e 

« 

contain absolute errors (Appendix C) which the ATS-6 model cannot account for. 

However, our previous work using landmark measurements has provided us with 
bounds on these errors and the rate at which they change as a function of 
image coordinates. The worst— case estimates correspond to about +2.5 pixels 
(+ 3.6 X 10~*^ radians). The errors tend to be oscillatory as a function of 
line or element position with a period of 100 to 200 pixels. Thus for 
typical feature displacement measurements (corresponding to tens of pixels), 
e can be considered to be a constant locally. (It would be appropriate to 
remind the reader at this point that the absolute mirror scan nonlinearities 
as a function of image coordinates repeat from image to image.) 

Taking the transpose of (19) and operating on (17) and (18) 


LVl 


— T'' 
= R r. 


PFl 


= (I + e)R,R r. 


1 1 LVl 


= r. 


LVl 


+ r 


(20a) 


^^LV2 ^ ^PF2 ^^^1^1^LV2 ^^LV2 ’ 


where 


r “ er. 
e 


LV * 


(20b) 


( 21 ) 


Application of R then, results in the correct LV vectors to within 
the same additive vector constant. It is shown in Appendix C that this 
constant has a* negligible effect in computing cloud velocities or displacements. 


ORIGINAL PAGE IS 
OF POOR QUALITY'/ 
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In addition location errors on the earth due to r^ is given in this 
appendix. 


4. General Case of n-image Sequence 

The case for images , T^, ... simply a generalization of the 

2- image case. The error-free LV to PF frame matrix for T, can be related 
to by the image pair relative-attitude transformation matrices 
where 

R = (I + E . ) ••• (I + E„^)(I + E,^)R, 

n n-l,n ' 23' 12 1 




(22) 


where only first order terms are retained. Thus the analog to (10), generalized 
to n images j is 

^n-1 


'r * = r 

PFn PF 




’^PFl \ ^LVn ’ 


where in practice is computed by evaluating the expression 


r =; T 

PFn PF 


fn - <23) 

' X=1 X=1 ' 


where L, .,,(L,,_), E. ,,.(L.,-) are evaluated from the T, to T - earth-edge 
X , x+1 x+1 X , 1+1 x+1 X x+x 

displacement measurements evaltmted at where Is the number 

of the "center of gravity" of the feature being tracked. 


Application of (19) to the r’ . results in 

IT r i 


IT* “ 3T IT 

LVi ^LVi e 


1 ^ ^ ^9 • • • H 


(24) 
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APPENDIX C. EFFECT OF ABSOLUTE ATTITUDE ERROR ON THE ACCURACY OF WIND 
MEASUREMENTS AND FEATURE LOCATION 


As discussed in the main body of the report and Appendix B, the 
transformation from the ATS— 6 image frame (PF) to the Local Vertical (LV) 
contains an absolute error term because of landmark measurement errors 
and absolute mirror scan nonlinearities. This means that the computed 
location in earth coordinates of a feature on the ATS— 6 image will be in 
error. Because this error is the same for successive images and changes 
slowly as a function of image position compared to features such as clouds 
"moving" on these successive images, the error of the measured displacements 
and consequently the errors in the wind measurements will be negligible 
but the errors in absolute location resulting from the absolute attitude 
uncertainty error will become more important. The purpose of this appendix 
is to compute the velocity and location errors as a function of 
angle measured from the SSP due to the absolute attitude error. 

From eqs. (4) and (20) in Appendix B, 






(tj) + 


r 

e 


( 1 ) 


“LV 


(tj 




H2) 


where 








-> 

r 

e 





= correct LV lanit vector which points from the satellite 
to a cloud at time t^ 

= computed LV unit pointing vector for cloud at time t^ 

= vector error resulting from absolute attitude error 
= computed LV unit pointing vector for cloud at t^ 
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“ change in computed pointing vectors due to clouds 
motion in the time interval At^^^ from t^ to t^ , 

where [!?J , |A?^^|t « |lrj^^(t^)[| = i. 

The geometrical relationship of these terms is shown (greatly exaggerated) 
in Figure C.l. where it is assumed all vectors are coplanar. Figure C.2. 
shows how these vectors relate to the satellite/earth system. 

In Figure C.2. the following are defined: 

P = point on earth’s surface of cloud at time t^ 

9 = angle at earth center between SSP-Satellite line and cloud 

(j) = angle between SSP-Satellite line and r .,(t ) 

j-»v 2 

R = earth radius = 6378 km 
K = 6.61 for geosynchronous satellite. 





FIGURE C.l. Geometrical 
Relationship Among Pointing 
Vectors 



FIGURE C.2. Satellite/Earth 
Geometry 
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Clearly the greatest error resulting from occurs when it lies in 
the plane formed by and r^^Ct^) and points in the same direction 

as Since an upper bound on the errors is desired all subsequent 

analysis will be confined to a plane and it will be assumed that 

point in such a direction as to increase <j>. With these stated conditions ' 
the computed velocity V can be treated as a scalar. 

The correct velocity V of the cloud from times t^^ to t^ at point P 
is given by 

V = C[6((i) +'A(j)) - 0((1))] , (3) 

where 


C = R/At 


12 


A(j> = change in 4> <iue to cloud motion 


6((j)) = value of 6 for a given value of <|) to be derived later on. 


Assuming A(J> « ^ and expanding 0 about keeping first-order terms, (3) 
becomes 


V = C 0’ A(j), 0' = d0/d^ . (4) 

The error in V, AV resulting from r^ is 

AV = C{[e(({) + A(^ + A^^) - 0(^ + A^)^)] - [0((j) + A<j>) - 0(^)]} 

= G[6’(<i> + AA ) - 0(<j>)] A<|. 

= C A(p A<t> , (5) 

where (5) was expanded as above keeping first-order terms and 




d^0~ 

de^J<t> 


A(p 


= \\r_ 


= error in resulting from vector error r 
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At this point it is convenient to compute 0', 9". From Fig. C.2. 
it is easily seen that 


tan<{» 


sin6 

K - COS0 


( 6 ) 


From (6) 0', 6" are computed and (6) is used to eliminate (jr, thus: 

, (7, 

e"(6) . e'(6) . 


( 8 ) 


where 0’ = 


d0 

d<f 


0 " 


d9* 

d4> 


The fractional error in velocity AV/V as a function of 0 can be expressed 
using (4), (5), (7) and (8) as 


V 


9"(9) 

e'(0) 


A* 


K sin0 (2 -f e* (9)) 
K COS0 - 1 


Acj), 


(9) 


Values of AV/V can be computed as a function of 0 when A^^ is assigned a 
value. As discussed before the worst-case estimate in the absolute attitude 
error is 

+ 2.5 pixels = + 2.5 pixels (1/120 ATS-6 pixels/deg) 

= + .021® = + 3. 6 X 10“^ radians. 

Using A<})^ =3.6 X 10“**^ radians in (9), we find that for 0 = 60“ , 

AV/V < .02 . 

Since AV/V increases with increasing angle 0 and recalling that 
almost all wind-velocity measurements are confined to study areas within 
an angle 0 of 50“ to 55“, we conclude that the effect on wind" velocity 
measurement errors due to absolute attitude- angle uncertainty errors is 
entirely negligible. 
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Location Errors 

The absolute location error S in earth coordinates (degrees and 
kilometers) is 

S(km) = R[0(0 + “ 0(^)1 

= R 0’(0) A<j)^ , (10) 

S (degrees) = 0'(0) A(}i . (11) 

e 

Using the values R = 6378 km and =3.6 X 10~^ and the expression for 0* (0) 
given by ( 7 ) 5 and substituting in (10) and (11) we arrive at the following table. 

Absolute Location Error 


Angle 
0 (degrees) 

for +2.5 Pixel Uncertainty 
+ AS (degrees) + AS (kilometers 

0 

0.12 

13.0 

10 

0.12 

13.3 

20 

0.13 

14.4 

30 

0.15 

16.3 

40 

0.18 

19.7 

50 

0.23 

25.8 

60 

0.34 

38.6 


TABLE C.l. Absolute Location Errors as a Function of 
Angular Distance from the SSP. These errors result from 
absolute landmark measurement and mirror scan nonlinearity 


errors . 
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APPENDIX D. PIXEL MEASUEEMENT ACCURACY REQUIRED FOR ACCURATE WIND MEASUREMENTS 


This appendix estimates the image-matching accuracy required between 
successive ATS-6 images to obtain wind-speed measurement accuracy to within 
+ 1 m/s. 

Referring to Fig. C.2. in Appendix C, the following are defined (it 
is assiamed that the image-matching error angles lie in the plane of the page) . 

A<j) — pointing angle error resulting from image-matching error 
between two successive ATS-6 images 

AP = pointing error in pixels corresponding to 

At = time difference between two -successive ATS-6 image (24.6 min) 


AV = velocity measurement error resulting from Acj) 

R = radius of earth (6378 km) 

K = 6.61 for a geosynchronous satellite. 

An image-matching error in A(j), results in a corresponding error in 
A6 and hence in the velocity V as follows : 

AV = R = R 7^ at the SSP, 


and 


AV 


d0 
At d<j) 


for an arbitrary 8, 


where 


( 1 ) 


AP (pixels) 7T_ rad 

-ATS-6 pixels , 180 deg 
degree ^ 

(K - cos6)^ + sin^9 
K COS0 - 1 


( 2 ) 


d<{) 


e 


eq. (7), Appendix C 


(3) 
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Substitution of (2) , (3) into (1) , solving for A? axid assuming the velocity 

\ 

error to be AV = + 1 meter/sec yields: 


AP = + (1,55 pixels) 


K cosQ - 1 


(4) 


(K — cos0)^ + sln^S 
for an ATS-6 two image sequence with a velocity error of +1 m/s. 


Equation (4) is used to generate the following table for a two- and three- 
image sequence. (The image-matching measurements between the first and 
second and the second and third images in a three-image sequence are assumed 
to be independent; thus the accuracy in a three-image sequence need only 
be l/VT times as accurate as for a two-image sequence to obtain wind-speed 
measurements accurate to 1 m/s.) 


Angle 
0 (degrees) 


Accuracy in Pixels AP Required to Obtain Wind 
Speed Accuracies Within + 1 m/s 

Two-Image Sequence Three-Image Sequence 

+ AP (pixels) + Vi~ AP (pixels) 


0 .29 .40 

10 .28 .39 

20 .26 ^ .38 

30 .23 .32 


40 


.18 


.26 


50 .14 .21 

60 .09 .13 


TABLE D.l. Accuracy in Pixels Required to Obtain Wind-Speed 
Accuracies to + 1 meter per second as a Function of Angle 
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APPENDIX E. EFFECT OF YAW ANGLE ERRORS ON WIND-MEASUREMENT ACCURACY 

The earth-edge shift measurement scheme used to account for relative 
attitude changes accounts for roll and pitch changes but cannot account 
for yaw-angle errors (a yaw-angle change results in an effective rotation 
of the earth image about the SSP) . This appendix derives expressions to 
estimate the errors on wind-velocity and wind-speed calculations due to 
yaw angle errors and computes wind-speed errors as a fmction of angle measured 
from the SSP line. 

For simplicity and without loss of generality, let the earth-centered 
XYZ coordinate system coincide with the earth-fixed axes relative to which 
latitude and longitude are measured (Fig. E.l.) and let the satellite be 
located on the + X axis. 



FIGURE E.l. Yaw-Angle Geometry 
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Define the following: 

r =» R[X,Y,Z] = vector locating a point on the earth, radius 

R, whose scaled components are X, Y, Z. "T" is the transpose 
of _the row vector. 

(9,^) = latitude, longitude corresponding to r 

A0y = yaw-angle change between two successive ATS-6 images 
separated in time by At 


(AU,AV) = apparent velocity errors in U- and V-components due to yaw- 
angle change A9^ 

0 = angle between X-axis and r 

AVy = apparent wind-speed error. 

The vector r expressed in 6 and <j) Is 

-> T " T 

r R[X,Y,Z] = R[cos0cos^, cos9sin<}», sin6] . (1) 


Now a yaw— angle change A0y between two consecutive image scans of the point 
located by r (in the first image) over a time At results in an effective 
TJ- and V-component velocity error given by the expressions: 


AU = R COS0 
At 


AV = R ^ , 
At ’ 


where A9,A(J) are the latitude and longitude errors resulting from A9y, 
The apparent change in r,Ar over the time interval At due to AS^ 
can be expressed in terms of an infinitesimal rotation by the angle A9 


( 2 ) 

(3) 

y 


about the X-axis: 
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A? 



1 

o 

o 

o 

! 


‘x' 

= R 

<D 

<3 

O 

O 


Y 


y 




0 -A0 0 

•- V 


-Z. 


y 


= R[0, A0 Z, -A9 Y]^ 

y y 


= R[0, A0^sxn0, — A0^cos9siii(i),] ' 


( 4 ) 


Now the relationship (1) can be differentiated with respect to 0,(j> and the 
differentials replaced by A’s to yield 

dr = Ar 


= R[-A({icos0sin(j) - A9sin9cosi}i, A(cos0sin<f>), A6cos9] 
The equality of the X-components in (4) and (5) yields 


T 


(5) 


A4COS0 = - 4J ^ n0 ^ s<j. 

sin<{» 


and the equality of the Z-components yields 


A9 = -A0ySin(j) . 


Substituting (7) into (6) results in 


( 6 ) 


(7) 


AiJ)cos9 = A9 sin0cos^ 


( 8 ) 


Substituting the EHS of (8) into (2) and the EHS of (7) into (3) results 
in the apparent U- and V-velocity component errors resulting from yaw-angle 
errors as a function of latitude and longitude; 


A0 

AU = R sin0cos(f) 

(9) 

A0 

AV = -R — sin<f. . 

(10) 
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The total wind-speed error AV is 

y 


A6 


AV = N/atJ^ 4- AV^ = R ^ \/sin^9cos^di + sin^i 
y At ^ ^ 


RA0 RA0 

- ... -. y . - cos^0cos^4> = - . ^ Vl - cos^0 . ' 

At At A 


RA0 

y 


At 


sinQ , 

A 


( 11 ) 


where is the azimuthal angle between the X-axis and r. ((H) could have 
been derived directly from very simple geometric considerations, but for the 
sake of completeness, it was desired to derive the AU and AV error terms.) 

Using the worst-case estimate of A9^ as + 2.5 pixels, the following can 
be evaluated: 


Aft = - 2.5 pixels IT rad 

y ^20 ATS-6 pixels ' 180 deg ’ 
degree 

R = 6.378 X 10^ m , 

At = 24.6 min (60 sec/min) , 

substitution of these values into (11) yields 

AVy = 1.57 sin0^ (m/s) , (12) 

for a 2- image ATS-6 sequence, and 

AV = sin0 =1.11 sin0. (m/s) (13) 

y vT 

for a 3-image sequence. From (12) and (13) Table E.l is prepared 



Angle 


Velocity Error AV due to + 2.5 pixel Yaw 
Error ^ 


e 


A 


(deg) 

2-image Sequence 
(+ m/ s) 

3-image Sequence 
(+ m/ s) 

0 

0 

0 

10 

.27 

.19 

20 

.54 

.38 

30 

.79 

.56 

40 

1.01 

.71 

50 

1.20 

.85 

60 

1.36 

.96 


TABLE D.l. Velocity Error as a Function of Yaw-Angle Error and Angle 


The low-level grid-point interpolated winds measured from ATS-6 images and 
coinciding with those winds measured from SMS-1 images (TABLE II. 2.) span 

an angular range about the ATS-6 SSP line from 22° to 40°. Thus 

the largest ATS-6 wind— speed error in a 3-image sequence that could result 

from a + 2.5 pixel yaw-angle error would range from .4 m/s to ,7 m/s. As 
seen from TABLES II. 1. and II. 2. these values are less than the rms wind- 
speed errors measured. Thus we can conclude for this image area studied 
that the effect of yaw-angle error on wind-speed measurements is small 
compared to other sources of error. 
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APPENDIX F, EFFECT ON EARTH-EDGE DISPLACEMENT MEASUREMENTS DDE TO 
ECCENTRICITY OF ATS-6 ORBIT 

As discussed in Appendices A and B, tlie effective attitude changes 
between successive ATS-6 images as a function of scan line is determined 
from the measured right and left earth-edge displacements between the 
successive image scan lines. 

Clearly for nonzero eccentricity of the satellite's orbit, the angular 
size of the earth as seen from the satellite will change from one image 
to the next. 

If this change were large enough, it could affect the attitude-change 
estimates between consecutive images derived from the earth-edge measurements. 
This appendix places an upper-bound estimate on the angular change of the 
earth's image between consecutive ATS-6 images. 


Satellite 



FIGURE F.l. Satellite/Orbit Geometry 



F-2 


Let: 

r = distance from the earth-center to the-satellite 
R = radius of earth 
V = true anomaly angle 

6 = half-angle that the earth subtends at the satellite 

6 = eccentricity of ATS-6 satellite’s orbit = 2 X 10~‘’ 

A0 - change in 6 over a time interval At (for two consecutive ATS-6 

images At = .41 hours) due to the eccentricity satellite’s orbit 

a = semi-major axis of orbit = KR, K = 6.61 geosynchronous 
satellite 


The equation of the orbit in terms of r, V and e is: 


r 


P 

1 + ecosV 


P 


a(l - e^) 
KR(1 - ^ 


and since e << 1, P = KR and (1) becomes 


r = KR(1 - ecosV) . 


( 1 ) 


( 2 ) 


Now the greatest change in r for a change in V occurs when 

dr/dV = KRe sinV has its maximum value which occurs when V = Vm = 90°; 

thus 


(dr/dV) = KRe . (3) 

max 

From Fig. F.I., sin6 = R/r; differentiating both sides and replacing the 
differentials with A’s yields: 



r^ (KR)2 K^R 


( 4 ) 



F-3 


The maximum change in A9 is computed from (4) as follows; 


A6 


Ar 


max max 


= — fe) AV 


(5) 


^ AV = I AV 

r2r ^ 


( 6 ) 




where (3) was used to replace (’^j with KRe. 

^ ' max 

Now AV = (15 (♦ 41 hr) = 6.15 deg for an ATS-6 2-image sequence, thus 

(6) becomes 


9X1 0“** 1 

^0 = — 6.15 deg = 1.86 X 10"^^ deg 

max 6.61 

= (1. 86 X lO"*" deg) (120 ATS-6 pixels/deg) 

= .02 ATS-6 pixels. 

Thus the worst-case estimate of the eccentricity results in an angular 
change between successive ATS-6 earth images which is entirely negligible. 
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APPENDIX C. METHOD OF OBTAINING MIRROR- SCM-OFFSET CORRECTION CURVES 

The nature of the mirror-scan nonlinearity problem has been discussed 
in previous reports. In order to correct the alternate scan offset caused 
by this effect, we need a table of offset values as a function of element number 
(i.e. a AE(E) function). For our initial efforts to produce a AE(E) function, 
we used the McIDAS cloud-tracking program to measure the displacement of 
a feature seen in the odd number scans of an image to its position in the 
even numbered scans of the same image. The method described here is a 
somewhat more automated scheme and does not require viewing the ATS-6 images. 

The. method has been applied to an IR image (74195 173134Z) and gives good 
agreement with the old method with less scatter of individual points about 
a pol 3 momial least-squares fit curve (FIG. G. 1.). 

In the automated method of computing the AE(E) function a correlation 
value is computed for the match between a small segment of an odd scan and 
a shift segment of an adjacent even scan. The amount of shift which gives 
the best correlation for that small line segment is taken to be the AE 
value for that scan line and element location. Values are computed over 
many scans and elements, then averaged over the scans. The result is a table 
of AE values as a fimction of E (element). The table values are then smoothed 
by using a least squares fit polynomial. In more detail the method is as 
follows: 

Let: 

e - element number (1 e _< 2400) 
s = scan number (1 s ^ 1200) 

p^(s,e) = pixel digital value at a given even-scan, element position 



G-2 


p^(s,e) = pixel digital value at a given odd-scan, element position 
5 = an element shift value 

t = a small number defining the length of the scan segment used 
for a correlation 

C(s,e,S) = a measure of the match between line segments from adjacent 
odd and even scans. 

Now define C by: 

e+t 

C(s,e,6) = 2 Ep (s,e’) - p (s + 1, e' + 5)]^ (s odd) 

e’=e-t ° ® 

The element shift, AE, for a given point on the image (s,e) is the value 
of S that gives a minimum value to C: 

C(s,e,AE(s,e)) ^ C(s,e,5) V6 

After an array of these values is generated, a weighted average over scans 
is taken: 

S w(s,e)AE(s,e) 

AE(e) = 

2 w(s,e) 
s 

The weighting factor is the range of brightness values from the odd scan 
used in computing G. 

w(s,e) = MAX(p^,Cs,e’ )) - MIN(p^(s,e’)) 
where e“t_5.e’_^e + t. 

The reason for using this weighting is simply that high contrast features 
should, in general, produce a correlation value of more significance than 
features of nearly uniform brightness. 
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The computer program to do this used the values; t = 7 and +4 6 ^ 12. 

Values of AE(s,e) were computed for every eight scans for the middle third 

of the image (s = 399, 407, 799) and for every 10 elements across most 

of the width of the image (e = 100, 110, ..., 2300). Using this line and 
element range some AE values will be computed for points off the earth. When 
off-earth and on-earth points are averaged together the weighting factor will 
give only a very small or zero weighting to the off-earth points. For some 
element values no on-earth points are encountered. Vfe found that these 
points could be determined from the table of AE(e) values. There is a 
discontinuity between the on-earth and off-earth values. 

The AE(e) function is then used to create a corrected image by shifting 
even scans. An improved method of correcting the alternate scan offset has 
also been developed. In the old method a fixed shift value was used across 
the width of a McIDAS image (672 pixels). Although this gave fairly good 
results, there could be a one or two element alignment error near the edges 
of the image. The new method, as before, uses the odd scans as a fixed 
reference. However, the required amount of shift is computed, or looked up 
in a table, for each pixel in the even scans. Thus the amotint of shift varies 
from one side of the image to the other. Pixels are dropped or doubled, as 
appropriate, between regions of different shift values. 



ELEMENT OFFSET 
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RLTERNRTE SCRN OFFSET. RTS-6 IR 74195 173134Z. 
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FIGURE G.l. Mirror-Scan Nonlinearity Curve AE(e) using newly developed procedure. 






